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Using the recently introduced model NEXUS 3, we calculate for pp,np 
and pp collisions the excitation function of particle yields and of average 
transverse momenta of different particle species as well as rapidity, xp 
and transverse momentum distributions. Our results are compared with 
available data in between y/s = 5 GeV and 65 GeV. We find for all observ- 
ables quite nice agreement with data what make this model to a useful 
tool to study particle production in elementary hadronic reactions. 

1 Introduction 

The recently published Q and Q data in pp collisions at 1 58 GeV by the 
NA49 collaboration have invalidated all customary quark-diquark string 
models [T] EI , which have been employed since long time to describe pp as 
well as AA collisions. In all those models, the valence quarks are taken 
as the end points of the strings, which provokes that due to the string 
topology \2\ more Q than Q are produced, in disagreement with the data. 
This effect has been verified in detailed calculations. 

In the new NEXUS 3 model, the observed particles are produced by 
two sources: a) strings which are formed by sea (anti) quarks and which 
are therefore symmetric with respect to the exchange of a particle and an 
antiparticle and b) excited remnants, which decay statistically. Whereas 
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Figure 1 : Multiple elementary interactions (Pomerons) in NEXUS . 

the string produces an equal number of Q and Cl, the remnant favors due 
to its finite baryon number the production of Q, in accordance with the 
data. 

In this paper, we calculate the results of the NEXUS 3 model for parti- 
cle yields, rapidity distributions, transverse momentum distributions and 
average transverse momenta. Where experimental data are available we 
compare with them. We will show that the prediction of this model, once 
the necessary parameters are fixed to describe the CERN SPS pp data at 
158 GeV and the excitation function of the multiplicity of charged par- 
ticles, are in good agreement with the experimental data down to beam 
energies as low as 10 GeV . A similar, however less complete, comparison 
between data and string fragmentation models has recently be performed 
by H. Weber et al. ©. They compare the URQMD and HSD model with 
data of nonstrange or single strange particles. 

2 NEXUS 3 

NEXUS 3 is a self-consistent multiple scattering approach to proton- 
proton and nucleus-nucleus scattering at high energies. The basic fea- 
ture is the fact that several elementary interactions, referred to as Pomerons, 
may happen in parallel. We use the language of Gribov-Regge theory to 
calculate probabilities of collision configurations (characterized by the 
number of Pomerons involved, and their energy) and the language of 
strings to treat particle production. 

We treat both aspects, probability calculations and particle produc- 
tion, in a consistent fashion: In both cases energy sharing is considered 
in a rigorous way |4|, and in both cases all Pomerons are identical. This 
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Figure 2: Inelastic scattering in pp. Partons from the projectile or the 
target proton interact via elementary interactions (the corresponding pro- 
duced particles being represented by horizontal lines) , leaving behind two 
remnants. 



is one new feature of our approach. Another new aspect is the necessity 
to introduce remnants: The spectators of each baryon form a remnant, 
see Fig. [TJ They will play an important role on particle production in 
the fragmentation region and at low energies [E Lab =40-200 GeV). In the 
following we discuss some more details of our approach. 

We first consider inelastic proton-proton scattering. We imagine an 
arbitrary number of elementary interactions to happen in parallel, where 
an interaction may be elastic or inelastic, see Fig. El The inelastic am- 
plitude is the sum of all such contributions in with at least one inelastic 
elementary interaction is involved. 

To calculate cross sections, we need to square the amplitude, which 
leads to many interference terms, as the one shown in Fig. E3a), which 
represents interference between the first and the second diagram of Fig. 
El We use the usual convention to plot an amplitude to the left, and the 
complex conjugate of an amplitude to the right of some imaginary "cut 
line" (dashed vertical line). The left part of the diagram is a cut elemen- 
tary diagram, conveniently plotted as a dashed line, see Fig. E3b). The 
amplitude squared is now the sum over many such terms represented by 
solid and dashed lines. 

When squaring the inelastic amplitude, all of the terms which corre- 
spond to the same final state interfere. For example, a single inelastic 
interaction does not interfere with a double inelastic interaction, whereas 
all the contributions with exactly one inelastic interaction interfere in- 
dependent of the number of elastic collisions. So considering a squared 
amplitude, one may group terms together representing the same final 
state. In our pictorial language, this means that all diagrams with one 
dashed line, representing the same final state, may be considered to form 
a class, characterized by m — 1 - one dashed line ( one cut Pomeron) - and 
the light cone momenta x + and x attached to the dashed line (defining 
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Figure 3: Inelastic scattering in pp. a) An interference term of cross 
section, b) Represented with our simplified notations. 




Figure 4: Class of terms corresponding to one inelastic interaction. 



energy and momentum of the Pomeron). In Fig. |4l we show several dia- 
grams belonging to this class, in Fig. El we show the diagrams belonging 
to the class of two inelastic interactions, characterized by m = 2 and four 
light-cone momenta 

Generalizing these considerations, we may group all contributions 
with m inelastic interactions. The sum of all these terms represents the 
probability of having m inelastic interactions with x\..x\ m , x{...x^ m at a 
given impact parameter. Integrating over impact parameter provides the 
corresponding cross section. By this we obtain a probability distribution 
for the number of elementary interactions (number of Pomerons) and the 
momenta of these Pomerons. 

How to form strings from Pomerons? No matter whether single-Pomeron 




Figure 5: Class of terms corresponding to two inelastic interactions. 
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Figure 6: Each Pomeron is identified with two strings. 




Figure 7: Remnants in single (two strings) and double scattering (four 
strings): in any case, two remnants contribute. 



or multiple -Pomeron exchange happens in a proton-proton scattering, all 
Pomerons are treated identically. Each Pomeron is identified with two 
strings, see Fig. |BJ 

The string ends are quarks and antiquarks from the sea. This dif- 
fers from traditional string models, where all the string ends are valence 
quarks. Due to the possibility of having a large number of Pomerons, 
this is impossible in our approach. The valence quarks stay in remnants. 
Being formed from see quarks, string ends from cut Pomerons have com- 
plete flavour symmetry and produce particles and antiparticles in equal 
amounts. 

Remnants are new objects, compared to other string models, see Fig. 
H The partonic content of a remnant is as follows: three valence quarks 
and the corresponding antiparticles of the partons representing the string 
ends. The masses of remnants are assumed to be small compared to the 
kinetic energies involved and are therefore neglected for the calculations 
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Figure 8: a) The most simple and frequent collision configuration has two rem- 
nants and only one cut Pomeron represented by two q — q strings, b) One of the 
q string-ends can be replaced by a qq string-end. c) With the same probability, 
one of the q string-ends can be replaced by a qq string-end. 



of multi-Pomeron configurations. To obtained finally the masses, one 
parameterize the mass distribution of a remnant as P(m 2 ) oc (m 2 )~ a , m 2 e 
( m min5 x+s )> where s is the squared energy at center mass system, 

^min IS 

the minimum mass of hadrons to be made from the remnant's quarks and 
antiquarks, and x + is the light-cone momentum fraction of the remnant 
which is determined in the collision configuration. Through fitting the 
data at 158 GeV we determine the parameter a = 1.5. Remnants decay 
into hadrons according to n-body phase space 

The leading order and therefore the most simple and most frequent 
collision configuration has two remnants and only one cut Pomeron rep- 
resented by two q — q strings as in Fig. |8^.. Besides the three valence 
quarks, each remnant has additionally a quarks and an antiquark to 
compensate the flavour. 

In NEXUS 3, this most simple approach is slightly modified by allow- 
ing with a small probability P qq that an antiquark q at one of the legs of 
the Pomeron is replaced by a diquark qq. The corresponding string ends 
are then a diquark and a quark. In this way we get quark-diquark (q — qq) 
strings from cut Pomerons. The qqq Pomeron end has to be compensated 
by the three corresponding antiquarks in the remnant, as in Fig. 
The (3q3q) remnant may decay into three mesons (3M) or a baryon and an 
anti-baryon (B + B), but the 3M mode is favored by phase space. For sym- 
metry reasons, the q leg of a cut Pomeron is replaced by an antidiquark 
qq with the same probability P qq . This yields a q — qq string and a (6q) 
remnant, as shown in Fig. |8k. The (6q) remnant decays into two baryons. 
Since q — qq strings and q — qq strings have the same probability to appear 
from cut Pomerons, baryons and antibaryons are produced in the string 
fragmentation with the same probability. However, from remnant decay, 
baryon production is favored due to the initial valence quarks. 
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With decreasing energy, the relative importance for the particle pro- 
duction of the strings decreases as compared to the remnants, because 
the energy of the string is lowered as well. If the mass of the string is 
lower than the cut off, it will be discarded. However, the fact that an 
interaction has taken place is taken into account by the excitation of the 
remnant which follows still the above mentioned distribution. 

3 Results 

To fix the main parameters for the Pomerons and the strings, we use 
the total cross section for pp scattering, the excitation functions for all 
identified charged particle yields from sqrt(s)=10 GeV to 2 TeV and the 
rapidity, pt and multiplicity distributions for at 200 GeV lab and cms and 
1800 GeV for charged particles. To adjust the parameters of the remnant 
decay and the multistrange baryon production, we mainly use the NA49 
results at 158 GeV lab. Then, with this set of parameters, we compare 
the simulation results to all pp data at quite wide energy range and also 
to np data at beam energy 40 GeV. This comparison work is not trivial, 
firstly because with decreasing beam energy the importance of the rem- 
nants for the particle production increases as compared to the pomerons. 
Secondly, pp collisions fix only the sum of the contributions of remnants 
and pomerons which contribute to different rapidity regions. The com- 
parison of the predictions for np collisions with data reveals whether each 
individual source of particle production is correctly described. 

3. 1 Hadron Multiplicities 

3.1.1 Energy Dependence of Average Multiplicities 

We start our investigation with the 4n multiplicities at different energies. 
In fig. |§]we display the excitation function of charged particles |7|-|25| as 
compared to the NEXUS 3 results. The A A and K s have been identified 
and their decay products are not included here. The agreement between 
calculation and experiment shows that this approach is able to describe 
excitation functions. 

Now we come to the excitation function of identified hadrons. The 
predictions for 4tc multiplicities of identified hadrons are shown in figs. 
110111 1 II 1 21 and here as well we find agreement with most of the data over a 
wide energy range, from y/s = 5 GeV to 63 GeV. The largest discrepancies 
we observe for the <p meson and for the p at low energies. It would be 
interesting to see in which rapidity region this discrepancies appears ( 
especially because the A is reasonable well reproduced but unfortunately 
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Figure 9: Excitation function of charged particle |7|-|25| as compared to 
the NEXUS results. The A A and K s have been identified and their decay 
products are not included. 

such data are not available. 

3.1.2 K/n ratios 

Special interest has recently gained the excitation function of the K/tt 
ratio because it has been claimed that the appearance of a maximum in 
the K + /n + excitation function, observed in heavy ion reactions, may be 
a signature that the plasma of quarks and gluons is formed. Although 
we are not concerned with heavy ion reactions here we display these 
ratios for later reference in fig. 1131 Neither data nor calculation show a 
maximum of the K + /n + ratio in this elementary collision. 

3.1.3 Rapidity Spectra 

We come now to more detailed information by studying the rapidity distri- 
bution. Fig. 1 141 shows the rapidity distribution of a multitude of hadrons 
for the reaction 40 GeV pp. Where data from the NA49 collaboration are 
available we have included them in the plot. The numbers give the av- 
erage multiplicity of the hadrons in 4ir. We see that the experimental 
data are reasonable well described. The non-strange baryons as well as 
those which contain one strange quark show a double hump structure, 
the others are peaked at midrapidity. This is a consequence of the three 
source structure (two remnants and Pomerons) in our approach. The 
leading baryon has still the quantum number of the incoming baryon but 
is moderately excited. Therefore it may disintegrate into baryons whose 
quantum numbers differ not too much. Also the calculation of the np re- 
action at 40 GeV, displayed in Fig. 1 151 reproduces the data quite nicely. 
We see of course a much lower number of protons in the fragmentation 
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Figure 10: Excitation function of the multiplicity of identified hadrons as 
compared with the available data (squares) |7|-|50|. 
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Figure 12: The same as Fig. \W[ 
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Figure 13: Theoretical and experimental K/n ratio. Squared points are 
from 1 45 1, triangle points represent the K s /ir ratio and are from |26|-|44|, 
compiled in Q. 



region. This difference to the pp collision is well reproduced. This agree- 
ment validates the correct description of the basic mechanism for par- 
ticle production in the fragmentation region. A comparison between the 
spectra shows that the difference between pp and np in the pion spectra 
extends to negative rapidities (as seen as well in the data). For A and A 
as well as for the other charge zero particles, the difference between pp 
and np is negligible. Naturally ) are more copiously produced if the 

projectile is a p(n). 

Fig. 1161 shows the rapidity spectra for pp at 100 GeV. Again we see a 
quite reasonable agreement between data and calculation. At 158 GeV, 
Fig. 1 1 71 we included the experimental A, A, E, E spectra which have 
recently been published 1 1 1 . We observe also more Q than Q as seen in 
experiment 1 6 1 . This is a consequence of the modification of NEXUS 3 
explained in (2 as compared to the original NEXUS 2 version |4| which 
yields more Cl than Q due to the string topology . 

For ISR energies there are only rapidity distributions for a given trans- 
verse momentum. In figs. 1181 and \W\ we display for p T = 0.4 GeV/c the 
rapidity distributions for n + , ir~, K + , K~, p, p for energies in between E cm 
= 23 GeV and 63 GeV. We see that also here the spectra agree well with 
the data where data are available. 

Fig. 1201 shows the longitudinal xf distributions of identified hadrons 
from pp collisions at a beam energy 400 GeV. In order to be comparable 
with the before-mentioned rapidity distributions, we use a logarithmic 
representation of the x-axis. We see that also here the spectra agree 
reasonably well with the LEBC-EHS Col. data|50|. Please note that there 



12 




-3 3 -3 3 -3 3 -3 3 

rapidity y rapidity y rapidity y rapidityp/fp 




-3 3 -3 3 -3.5 3.5 -3.5 3.5 

rapidity y rapidity y rapidity y rapidity y 




-3 3 -3 3 -3 3 -3 3 



rapidity y rapidity y rapidity y rapidity y 




-3 3 -3 3 -3 3 -3 3 



rapidity y rapidity y 5 rapidity y g rapidity y 




-3 3 -3 3 -3 3 -3 3 



rapidity y rapidity y rapidity y rapidity y 

Figure 14: Rapidity spectra (lines) and A-k multiplicities (numbers) of iden- 
tified hadrons from proton-proton at beam energy 40 GeV. Data are from 
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Figure 15: Rapidity distribution (lines) and 4tc multiplicities (numbers) 
of identified hadrons for neutron-proton collisions at a beam energy of 40 
GeV. Data are from 1 6 1 . i a 
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Figure 16: Rapidity distribution (lines) and 4tc multiplicities (numbers) 
of identified hadrons for proton-proton collisions at a beam energy of 100 
GeV. Data are from |6|. i e 
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Figure 17: Rapidity distribution (lines) and 4n multiplicities (numbers) of 
identified hadrons for proton-proton collisions at a beam energy of 158 
GeV. Data are from 1 5) . i o 
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Figure 18: 7r + , tx and K + rapidity distributions at pt=0.4GeV/c for 
proton-proton collisions at ISR energies. Data are from 
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Figure 19: , proton and antiproton rapidity distributions at 

p t =0.4GeV/c for proton-proton collisions at ISR energies. Data are from 
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are two curves for the p's. This is due to the inconsistency of the two 
available data sets: z E dn/dxF and dn/dxF spectra. So dn/dxF spectra 
from our simulations are additionally plotted as dashed lines there. 

3.2 Transverse Momentum 

3.2.1 Average Transverse Momenta 

Fig. 12 II displays the excitation function of the average transverse momen- 
tum for all charged particles and several particle species in comparison 
with the experimental data|51 1. We see that with the exception of K + , the 
average transverse momentum is reproduced over the whole kinematical 
range. For the K + we underpredict the average transverse momentum at 
smaller beam energies. 

3.2.2 Transverse Momentum Distribution 

At SPS energies not only the mean transverse momenta but also the whole 
transverse momentum spectra is available and we compare these data 
with the NEXUS 3 predictions in figs. I22II23II24L Again, the agreement is 
quite reasonable up to p t = 2GeV. The calculation of the spectra at higher 
transverse momenta at SPS energies is beyond the limits of present day 
computers, because of the very small cross section of hard process. 

4 Conclusions 

Employing the recently developed NEXUS 3 model where the parameters 
have been mainly fixed to the pp data at E lab =158 GeV, we compare 
its predictions on average quantities, longitudinal and transverse spectra 
between y/s = 5 GeV and 65 GeV with the existing data of pp, np and pp 
collisions. 

We find a very nice agreement with data for A-k multiplicities, rapidity 
and x F distributions as well as for transverse momenta. This suggests 
that the basic mechanism of particle production is well described in the 
NEXUS 3 model. Based on this observation we can use this model to 
study pA and AA collisions which will be the subject of a forthcoming 
publication. 
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Figure 20: x F distributions of identified hadrons from pp collisions at 
beam energy 400 GeV. Data are from |50|. The vertical axis of dashed 
lines and the points in rho meson plots is dn/dx F . 
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Figure 21: Average transverse momentum of charged hadrons, pions, 
kaons and antiproton from 4n phase space at different energies. Data 
are from 1511 . 
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Figure 22: Transverse momentum spectra of proton, antiproton, charged 
pions and Kaons from pp collisions at E\ ah = 200GeV. Data are from 
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Figure 23: Transverse momentum spectra of proton, antiproton, charged 
pions and Kaons from pp collisions at Eu h = 300GeV. Data are from |52|. 
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pions and Kaons from pp collisions at Ei ab = 400GeV. Data are from |52|. 
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